/fb  A097  SO? 


echnical 

library 


TECHNICAL  REPORT  ARLCB-TR-79025 


THERMO-ELASTIC-PLASTIC  STRESSES  IN  HOLLOW  CYLINDERS 

DUE  TO  QUENCHING 

J.  D.  Vasilakis 
P.  C.  T.  Chen 


October  1979 

US  ARMY  ARMAMENT  RESEARCH  AND  DEVELOPMENT  COMMAND 

LARGE  CALIBER  WEAPON  SYSTEMS  LABORATORY 
BENffT  WEAPONS  LABORATORY 
WATERVLIET ,  N.  Y.  12189 


] 


AMCMS  No.  6111.01.91A0,0 


DA  Project  No.  1L161101A91A 
PRON  No.  1A-9-2ZA01-Y 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED 


DISCLAIMER 


The  findings  in  this  report  are  not  to  be  construed  as  an  official 
Department  of  the  Army  position  unless  so  designated  by  other  author¬ 
ized  documents. 

The  use  of  trade  name(s)  and/or  manufacturer (s)  does  not  consti¬ 
tute  an  official  indorsement  or  approval. 


DISPOSITION 

Destroy  this  report  when  it  is  no  longer  needed.  Do  not  return  it 


to  the  originator. 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Date  Entered) 

REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1  REPORT  NUMBER 

ARLCB-TR- 79025 

2.  GOVT  ACCESSION  NO. 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (and  Subtitle) 

S.  TYPE  OF  REPORT  &  PERIOD  COVERED 

THERMO- ELASTIC- PLASTIC  STRESSES  IN  HOLLOW 

CYLINDERS  DUE  TO  QUENCHING 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTHORfeJ 

J.  D.  Vasilakis 

P.  C.  T.  Chen 

8.  CONTRACT  OR  GRANT  NUMBERfe) 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Benet  Weapons  Laboratory 

Watervliet  Arsenal,  Watervliet,  NY  12189 
DRDAR-LCB-TL 

'  10.  PROGRAM  ELEMENT.  PROJECT,  TASK 

AREA  &  WORK  UNIT  NUMBERS 

AMCMS  No.  6111.01 .91A0 .0 

DA  Project  No.  1L161101A91A 
PRON  No.  1A-9-2ZA01 -Y 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

US  Army  Armament  Research  and  Development  Commanc 

12.  REPORT  DATE 

October  1979 

Large  Caliber  Weapon  Systems  Laboratory 

Dover,  New  Jersey  07801 

13  NUMBER  OF  PAGES 

-  20  

14.  MONITORING  AGENCY  NAME  ft  ADDRESS  (if  different  from  Controlling  Office) 

IS.  SECURITY  CLASS,  (of  thle  report ) 

UNCLASSIFIED 

is©  DECLASSI  F|  CATION/ DOWN  GRADING 

SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  thle  Report) 

Approved  for  public  release;  distribution  unlimited. 

,  17.  DISTRIBUTION  STATEMENT  (of  the  ebatrect  entered  In  Block  20,  If  different  from  Report) 

18.  SUPPLEMENTARY  NOTES 

Presented  at  the  25th  (Silver  Jubilee)  Conf  of  Army  Mathematicians, 

6-8  Jun  79,  Johns  Hopkins  Univ,  Baltimore,  MD. 

To  be  published  in  the  proceedings  of  the  25th  Conf  of  Army  Mathematicians. 

19  KEY  WORDS  (Continue  on  reveree  aide  If  neceeaary  and  Identify  by  block  number) 

Finite  Difference  Method 

Transient  Temperatures 
Thermo-Elastic  Plastic  Stresses 
Transformation  Stresses 

20.  ABSTRACT  (XZontbxue  am  revere*  mfdm  If  ne^eaeary  and  Identify  by  block  number) 

A  hollow  cylindrical  tube  rapidly  quenched  for  the  purpose  of  developing 
a  high  strength  material  structure  is  analyzed.  The  quenching  creates  severe 
thermal  stresses  early  in  the  quenching  cycle  while  later  the  material  trans¬ 
formation  by  virtue  of  a  volume  change  in  the  transformed  material  causes 
large  transformation  stresses.  The  transient  temperature  distributions  and 
the  elastic  treatment  of  the  stresses  has  been  treated  previously.  The  pres¬ 
ent  work  is  an  attempt  to  consider  the  thermo-elastic-plastic  aspects  of  the 

DD  i  jam *73  M73  EDITION  OF  t  NOV  65  IS  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dete  Entered ) 


S E CU  R1 TY  CLASSIFICATION  OF  THIS  PAG E(When  Data  Entmrmd) 


Cont  from  Block  20 

problem.  The  von  Mises!  yield  criterion  and  the  Prandtl-Reuss  stress  strain 
relations  are  used.  Results  are  calculated  based  on  a  new  finite  difference 
approach . 


SECURITY  CLASSIFICATION  OF  THIS  PAGEO*7ien  Data  Entered) 


TABLE  OF  CONTENTS 

Page 

SUMMARY  i 

INTRODUCTION  1 

PROBLEM  DESCRIPTION  2 

THERMAL  EQUATIONS  3 

STRESS  EQUATIONS  4 

NUMERICAL  COMPUTATIONS  6 

RESULTS  AND  DISCUSSION  8 

REFERENCES  10 

ILLUSTRATIONS 

1.  Residual  Stresses  in  a  Solid  Cylinder  Due  to  Material  11 

Transformation  (Transformation  Beginning  on  Outside 
Diameter  and  Progressing  Toward  Center) . 

2.  Residual  Stresses  in  a  Hollow  Cylinder  Due  to  Trans-  12 

formation  (Transformation  Occurring  Symmetrically 

from  the  Inner  and  Outer  Diameters) . 

3.  Residual  Stresses  in  a  Solid  Cylinder  Due  to  13 

Quenching . 

4.  Residual  Stresses  in  a  Hollow  Cylinder  Due  to  14 

Quenching. 


i 


SUMMARY .  A  hollow  cylindrical  tube,  rapidly  quenched  for  the 
purpose  of  developing  a  high  strength  material  structure,  is  analyzed. 
The  quenching  creates  severe  thermal  stresses  early  in  the  quenching 
cycle  while  later  the  material  transformation  by  virtue  of  a  volume 
change  in  the  transformed  material,  causes  large  transformation 
stresses.  The  transient  temperature  distributions  and  the  elastic 
treatment  of  the  stresses  has  been  treated  previously.  The  present 
work  is  an  attempt  to  consider  the  thermo-elastic-plastic  aspects  of 
the  problem.  The  von  Mises  yield  criterion  and  the  Prandtl-Reuss 
stress  strain  relations  are  used.  Results  are  calculated  based  on  a 
new  finite  difference  approach. 

I.  INTRODUCTION.  Watervliet  Arsenal  has  recently  been  develop¬ 
ing  techniques  for  the  production  of  large  caliber  weapons  using  a 
rotary  forge.  Force  hammers,  evenly  spaced  at  90°  intervals,  strike 
the  outside  diameter  of  a  hot  (1500°F-1600°F) ,  hollow  cylindrical 
preform  at  the  rate  of  200  blows/minute.  The  final  outside  tube 
profile  is  programmed  into  the  forge  itself,  and  the  wall  thickness  of 
the  tube  is  varied  as  preprogrammed.  The  inside  diameter  of  the  tube 
is  maintained  constant  by  a  mandrel  which  is  water  cooled.  After  the 
tube  has  been  formed,  it  is  allowed  to  cool  to  room  temperature. 

Once  formed,  the  tube  must  then  be  heat  treated.  This  procedure 
begins  by  heating  the  tube  to  1650°F  in  an  austenitizing  furnace  so 
that  the  austenite  phase  is  developed  throughout  the  material.  The 
tube  is  then  rapidly  quenched  so  that  the  desired  martensite  phase 
is  developed.  This  quenching  is  accomplished  by  spraying  a  large 
volume  of  water  on  both  the  inside  and  outside  diameters.  The  tube 
is  finally  put  through  a  tempering  furnace  at  1200°F. 

Interest  in  the  analytical  studies  of  the  process  first  arose  when 
cracks  began  developing  in  the  tube  during  quenching.  While  the  possi¬ 
ble  causes  of  quench-cracking  are  many,  most  often  they  are  associated 
with  the  material  used.  It  was  also  decided,  however,  to  look  into  the 
transient  temperatures  during  the  quenching  process  and  the  thermal  and 
transformation  stresses  involved.  The  transformation  stresses  occur 
mainly  due  to  volume  changes  in  the  material  as  it  transforms  from  one 
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phase  to  another.  As  they  are  due  to  volume  changes,  transformation 
stresses  can  be  treated  in  a  manner  similar  to  the  thermal  stresses. 
Although  quench- cracking  cannot  be  predicted  from  a  study  such  as  this, 
a  better  understanding  of  the  quenching  procedure  will  emerge  and  the 
relative  severity  of  different  quenching  procedures  would  be  known. 

The  transient  temperatures  and  the  zones  of  transformed  material 
assuming  a  linear  relationship  for  the  change  in  volume  between  the 
martensite  start  and  finish  temperatures  were  treated  in  [1] .  This 
reference  also  considers  the  thermal  and  transformation  stresses  assuming 
the  stresses  remain  elastic.  The  present  work  seeks  to  incorporate  an 
elastic-plastic  stress  analysis  into  the  problem.  In  view  of  the  previous 
results,  this  assumption  is  more  realistic*  The  temperature  and  stress 
problem  are  considered  uncoupled. 

II.  PROBLEM  DESCRIPTION.  The  problem  being  considered  is  that  of 
the  elastic-plastic  stresses  developed  during  the  quenching  process. 

These  stresses  are  due  to  both  the  transient  temperatures  that  exist  and 
the  transformation  stresses,. 

Most  of  the  elastic-plastic  analysis  work  on  thick-wall  cylinders 
concerns  itself  with  mechanical  loadings,.  Bland  [2]  does  consider  ther¬ 
mal  loads  on  a  thick  wall  tube0  Tresca's  yield  criterion  and  its  associ¬ 
ated  flow  rule  wereused  to  obtain  solutions  to  tubes  of  work-hardening 
material  subjected  to  both  internal  and  external  pressures.  The  temper¬ 
atures,  however,  are  steady  state, and  the  thermal  stresses  due  to  this 
steady  state  temperature  distribution  are  first  calculated  and  assumed 
elastic.  External  or  internal  pressures  are  then  applied  until  some 
desired  plastic  state  is  arrived  at.  S.  C.  Chu  [3]  used  the  incremental 
approach  for  solving  the  problem  of  elastic-plastic  thick-walled  tubes 
subject  to  transient  thermal  loadings.  The  von  Mises  yield  criterion  and 
Prandtl-Reuss  equations  are  used0 

In  the  area  of  elastic-plastic  analysis  for  transformation  stresses, 
the  bulk  of  the  work  comes  from  a  series  of  papers  by  Zwicky,  Landau, 
Weiner,  and  Huddleston  [4-6] .  Of  those  that  consider  the  cylinder  config¬ 
uration,  Weiner  and  Huddleston  [5]  used  the  Tresca  yield  criterion  and  the 
associated  flow  rule  to  compute  the  residual  stresses  in  the  cylinders. 

The  problem  for  the  transformation  stresses  was  solved  by  assuming  that 
the  volume  expansion  of  the  transformed  material  was  equivalent  to  that 
of  a  temperature  discontinuity  progressing  inward  from  the  surface.  They 
considered  a  solid  cylinder  of  incompressible  material.  Landau  and 
Zwicky  [6]  solved  a  similar  problem  using  the  von  Mises  yield  criterion 
and  its  associated  flow  rule.  They  assumed  a  compressible  material,  the 
yield  point  stress  to  be  a  function  of  temperature,  and  included  the  com¬ 
putation  of  transient  thermal  stresses. 
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The  problem  considered  here  is  that  of  determining  the  thermo- 
elastic-plastic  stresses  and  transformation  stresses  in  a  cylinder  due 
to  quenching.  The  thermal  program  developed  in  [1]  was  coupled  to  a 
Program  [7]  for  the  computation  of  elastic-plastic  stresses  in  a  thick- 
walled  cylinder  subjected  to  internal  and  external  pressure.  The  prob¬ 
lem  is  assumed  to  be  axisymmetric . 

The  computer  program  for  the  temperature  distribution  allows  for  a 
transient  analysis  with  temperature  dependent  material  properties  using 
an  implicit  finite  difference  scheme.  The  computer  program  for  the 
elastic-plastic  stresses  uses  an  incremental  approach.  It  has  been 
altered  to  include  stresses  due  to  thermal  loads.  The  von  Mises  yield 
criterion  is  used  with  the  associated  Prandtl-Reuss  flow  rule.  The 
material  is  assumed  compressible  and  is  capable  of.  work-hardening  although 
for  this  work  the  material  was  assumed  to  be  elastic-perfectly  plastic. 

-HI*  THERMAL  EQUATIONS.  The  partial  differential  equation  for  the 
temperature  (T)  in  a  thick-wall  cylinder  with  inner  radius,  a,  and  outer 
radius,  b,  is  given  in  dimensionless  form  by 

7^  [ktT>r  Hi  -  cmpd)  f-  a) 

where  r  is  dimensionless  radial  distance,  k(T),  c(T),  p(T)  are  dimension¬ 
less  thermal  conductivity,  specific  and  density,  respectively,  and  t  is 
dimensionless  time.  The  dimensionless  quantities  are  defined  as 


r 


T  = 


T-Tp 

T.-T 
1  o 


t 


(2) 


k(T)  =  k0K(T),  c(T)  =  c0C(T),  p(T)  =  PqR(T) 

and  r  is  the  radius,  T  is  the  temperature,  k0,  c0,  p0  are  thermal  conduc¬ 
tivity,  specific  heat  and  density  at  reference  ambient  temperature  T0 ; 

T^  is  initial  temperature;  and  t  is  time. 


The  boundary  conditions  are  written  as 


and 


9T 

3r 


h]T  =  "gl  at  r  =  a/b 


8T  .  _ 

87  '  h2T  =  "g2 


at  r=l  . 
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With  the  boundary  conditions  expressed  in  this  manner,  different  condi¬ 
tions  at  the  boundary  can  be  specified.  If,  e.g.,  gj  =  0  and  h1  ±  0 
and  finite,  then  a  convection  type  boundary  condition  exists  on  the  inner 
surface.  If  h i  was  very  large  and  g1  =  0,  then  T  =  0  is  specified.  If 
^2  an<^  2 2  are  Doth  large  and  not  equal,  then  the  temperature  T  =  §2^2 
is  specified  in  the  outer  surface. 


IV.  STRESS  EQUATIONS.  The  use  of  finite  difference  equations  to 
solve  the  thermo-elastic-plastic  stress  problem  requires  expressing  the 
equilibrium  equation  and  the  equation  of  compatibility  at  each  node  at 
which  the  finite  difference  equations  are  desired.  The  Prandt 1-Reuss 
flow  rule  is  used  to  eliminate  the  incremental  stresses  so  that  what 
results  is  a  matrix  for  evaluating  the  incremental  radial  and  tangential 
strains  at  each  node.  The  required  equations  follow,  written  in  dimen¬ 
sionless  form.  The  problem  is  treated  as  plane  strain. 


The  equation  of  equilibrium  is  written 


90j»  o-p— oq 

8r  r 


(3) 


where 


cr  (=  — )  is  the  dimensionless  radial  stress 


ae 

Oq (=  — )  is  the  dimensionless  tangential  stress 
°o 

and  aQ  is  the  yield  stress  in  tension,  and  the  compatibility  equation 


3r 


0  +  e9-er  =  q 


(4) 


where 


:e 


£q(=  E  — — )  is  dimensionless  tangential  strain 


£r(=  E  — )  is  dimensionless  radial  strain 


and  E/a0  is  yield  strain  in  tension  when  E  is  Young's  Modulus, 
compressibility  of  the  material  is  expressed  by 


The 


where 


£  =  aT  +  * 

£  =  —  (£r+£g)  is  mean  strain 


(5) 


(a  +afi+a  )  is  mean  stress 


3  '•■'r  0  “z 


K 


K(=  — )  is  dimensionless  bulk  modulus 
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a(  aTi)  is  dimensionless  coefficient  of  thermal  expansion 

and 

ez  =  0  for  plane  strain. 

Traction  free  boundary  conditions  are  used 

ar  =  0  at  r  =  a/b  and  r  =  1.  (6) 

It  was  desirable  to  write  the  finite  difference  equations  in  terms 
of  strain  alone,  hence,  the  stresses  in  the  equations  of  equilibrium  had 
to  be  expressed  in  terms  of  the  strains„  This  was  accomplished  by  mod¬ 
ifying  a  plastic  stress-strain  matrix  [8]  which  was  derived  by  inverting 
the  Prandtl-Reuss  equations.  The  inverted  Prandtl-Reuss  equation  is 


{da} 


[Dp]{de}  - 


EadT 

d-2v)a0 


{1} 


(7) 


where  the  stress  vector  is  {da}  =  {dar,  dCQ ,  daz}  ,  the  strain  vector 
{de}  =  {der,  de0,  0},  and  {1}  represents  a  unit  vector.  The  plastic 
stress-strain  matrix  [DP]  is  given  by 


[DP]  = 


pl-v  °r ' 2 
l-2v  "  S 
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1-2V 


V 

l-2v 


'a  ' 
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S 


1  -V 

l-2v 


(8) 


The  primed  stresses  are  deviatoric  stresses. 


a. ' 
i 


i  =  r,0,z. 


(9) 


At  each  node  during  a  computation,  the  von  Mises  yield  criterion 

1  [(Vae)2  +  (a0_az)2  +  (Var)2]  =  1  (10) 

is  checked  to  see  if  plastic  deformation  has  progressed  to  that  node.  If 
not,  the  stresses  remain  elastic  and  can  still  be  computed  using  (8)  by 
setting  the  deviatoric  stresses  equal  to  zero.  The  matrix  [DP]  then 
becomes  the  same  matrix  as  would  exist  if  linear  elastic  behavior  had 
been  assumed.  The  quantity  S  is  given  by 
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where 


s  -l^c1  +  %) 


(ii) 


a=  7aij'aij 
is  the  equivalent  stress  and 
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7iur 


OV* 


l  2 


ae ' *  <V2) 


H'  = 


da 

dEn 


(12) 


(13) 


is  the  slope  of  the  equivalent  stress/equivalent  plastic  strain  curve  and 
is  a  measure  of  hardening.  The  increment  in  equivalent  plastic  strain  is 
given  by 


2  P  P 

deP  *  3  dEijdEij  ' 


(14) 


Vo  NUMERICAL  COMPUTATIONS,  The  Crank-Nicolson  representation  for 
finite  differences  of  the  partial  differential  equation  governing  the 
temperatures  in  time  is  [1] 


[(a+iAr)ki+%>n+%]Ti+1>n+1  + 


+  [-(a+iAr)k.  .  .  -  (a+(i-l)Ar)k.  ,  .  -c.  ,p.  f. 

i+%,n+%  x-h,i\+h  i,n+Vi,n+% 


2Ar 


A  •)  (a+(i-%)Ar)]T. 

At  i,n+l 


+  [ (a+ (i-1) Ar)k.  j  .  ]T.  .  =  [-(a+iAr)k.  ,  .  ]T.  , 

i-%,n+V  i-I, n+1  1  v  J  i+%,n+V  i+l,n 


7  At  ^ 

*  [ta+iAr)ki.%,I,+^(a*(i-l)*r)ki.H>nA-ci>nt%pi_nt!st— )(a*(i-%)4r)]Ti  n 

*  t-(a*ti-1>Ar)ki-%,n*%lTi-l,n  •  (15> 

The  equation  is  solved  twice : 


1.  At  n+%  step,  allowing  k,p,c  etc.  to  take  on  the  values  at  t=n 

step. 

2.  The  new  temperatures  are  then  used  to  evaluate  k,c,p,  at  n+h 
step  and  the  set  of  equations  re-evaluated  for  the  temperatures  at  the 
n+1  step. 

The  computed  temperature  distributions  at  each  full  time  step  are 
saved  on  disk  and  eventually  called  in  when  required  by  the  stress  program. 
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The  finite  difference  equations  are  (for  solid  cylinder). 
Compatibility: 


'riA£e.  ,  *  (2ri-ri-l)A£9.  -  (ri-ri_i)ie:r 

i  i 

'ri(Eei-Eei  >  -  (ri-ri-i)'e0.-er.) 

1  1-1  i  i 


Equilibrium: 


(16) 


■r .  Aa 
i 


r.  ,  -  tri"ri-l)Aae.  *  (2rrri-l>AV  * 

1“1  1  1 

-r-  (a  -a  )-  (r.-r.  ,)( a  -aQ  1 
i  r .  r .  /  v  i  i-l'  ^  r  0  J 
ii-l  ii 

Substituting  the  Prnndtl-Reuss  equations  into  that  of  equilibrium 

-r  D(r,9)Ae  -r  D(r,r)Ae  +  [-(r  -r  )D(0,9)  +  (2r  -r.  jD(r,9)]AeQ 
i-1  i-l  1-1  1  1-1 

+  [-(r.-r  )D(0,r)+(2r.-r.  )D(r,r)]Ae 

i  i-l  i  i-l  r. 

i 

r  P(Y 

ri[0r.  -°r.]*tri-ri-lH°9.-0r.’*ri  tATi-ATi-ll 

1-1  1  11 

at  i  =  1  (zero  radius  for  solid  cylinder) 

-AeQ  +  Ae  =  eQ  -  e 
91  rl  91  rl 


(17) 


(18) 


(19) 


at  i  =  n  (or  outside  boundary)  ar  =  0  or 

D(r,0)Ae0  +  D(r,r)Ae 

n 


EaATn 


rn  1-2V 


(20) 


For  the  hollow  cylinder,  a  boundary  condition  similar  to  i  =  n  can  be 
written  for  i  =  1. 

The  solution  procedure  for  the  transient  temperature  problem  is  as 
follows.  The  temperature  problem  is  solved,  and  the  temperature  distri¬ 
butions  at  their  computation  times  are  stored  on  disk.  These  distri¬ 
butions  are  called  into  the  thermo-elastic-plastic  stress  program  one 
at  a  time.  The  corresponding  thermal  stresses  are  calculated  and  each 
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node  checked  to  see  if  the  yield  criterion  is  satisfied.  If  not,  the 
problem  is  still  assumed  to  be  elastic,  a  new  temperature  distribution 
is  called  in,  and  new  stress  increments  calculated.  The  stresses  are 
updated,  and  the  yield  criterion  checked  again.  When  the  stresses  at 
a  point  are  found  to  satisfy  the  yield  criterion  the  node  is  identi¬ 
fied,  and  the  stress  increments  at  that  node  from  the  next  set  of 
temperatures  are  computed  using  the  Prandtl-Reuss  equation  or  [Dp] 
matrix  identified  earlier.  This  procedure  is  continued  with  new  sets 
of  temperature  called  in  and  with  the  tracking  of  the  elastic-plastic 
boundary (s)  with  time.  The  resultant  stresses  that  exist  after  a 
steady-state  or  uniform  temperature  distribution  is  reached  are  the 
residual  stresses. 

The  solution  procedure  for  the  transformation  stresses  is  similar 
and  will  be  described  in  the  next  section. 

VI.  RESULTS  AND  DISCUSSION.  Several  runs  were  made  for  the 
stresses  due  to  the  transient  temperatures  and  for  the  transformation 
in  both  solid  and  hollow  cylinders.  For  the  results  presented  here, 
the  following  data  were  used: 

E  =  30xl06  psi,  aQ  =  30xl03  psi 

a  =  7 . 75x10-6/°F,  T^  -  1250°F  a  aT.  =  .0097 

v  =  .3 

h2  =  12.2,  hj  =  12.2  and  6.1 

The  first  results  shown  are  those  for  the  transformation  stresses. 
These  stresses  can  be  computed  using  the  thermal  stress  formulation  if 
one  replaces  aT  or  the  thermal  expansion  by  the  linear  expansion  of  the 
transformation.  If  the  material  expansion  due  to  the  transformation  is 
isotropic,  this  linear  change  is  1.3  the  volume  expansion.  As  an  exam¬ 
ple,  the  volume  expansion  in  going  from  the  austenite  to  the  martensite 
structure  for  steel  is  about  3%-4%„  In  the  quenching  of  a  solid  cylin¬ 
der,  the  transformation  begins  on  the  outer  surface  and  progresses 
inward  to  the  center.  Figure  1  shows  the  residual  stresses  in  a  solid 
cylinder  due  to  a  transformation  occurring  in  the  material.  The  insert 
shows  the  temperature  function  as  it  progresses  inward.  It  is  of  unit 
height  in  the  transformed  material  and  zero  in  the  untransformed  mate¬ 
rial.  The  transformation  is  assumed  to  be  occurring  over  eight  nodes 
or  8-s  of  the  cylinder  (indicated  by  N  in  Figure)  ,  and  a  linear  relation 
is  assumed  over  this  length.  Initially,  Oq  is  compressive  near  the 
outside  radius  when  the  transformation  just  begins  as  the  material 
wants  to  expand  but  is  prevented  from  doing  so  by  the  surrounding 
untransformed  material.  As  the  transformation  progresses,  however,  Oa 
slowly  changes  sign  and  becomes  tensile.  The  computer  run  was  stopped 
just  before  the  transformation  was  complete,  and  that  is  the  reason  for 
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the  behavior  of  the  stresses  near  the  bore.  The  transformation  at  r  = 

0  had  just  started  when  the  run  was  stopped. 

Figure  2  shows  similar  results  for  the  hollow  cylinder.  The 
assumption  is  made  that  the  transformation  progressed  evenly  from  the 
inside  and  outside  surfaces.  Tangential  stresses  on  the  outside  sur¬ 
face  again  were  initially  compressive  and  slowly  changed  to  tensile 
stresses  while  those  in  the  inside  radius  always  remained  compressive. 

Figures  3  and  4  show  the  residual  stresses  that  exist  due  to  the 
transient  temperatures  from  the  quenching  process.  Figure  3  represents 
the  results  for  the  solid  cylinder.  The  large  axial  stress  due  to  the 
plane  strain  assumption  is  easily  seen.  As  the  quenching  begins,  the 
outside  surface  cools  and  wants  to  contract.  It  is  prevented  from 
doing  so  by  the  surrounding  material  and  therefore  Oq  is  initially  a 
tensile  stress.  The  elastic-plastic  boundary  begins  on  the  outside 
surface  of  the  cylinder  and  moves  toward  the  center. 

Figure  4  shows  similar  results  for  the  hollow  cylinder.  The  fig¬ 
ure  shows  the  residual  stresses  when  equal  convection  type  boundary 
conditions  are  used  on  both  the  inside  and  outside  diameters.  These 
are  shown  by  the  solid  lines.  A  comparison  is  made  with  the  same 
problem  when  the  convection  boundary  condition  on  the  inside  diameter 
is  decreased  by  50%.  A  dotted  line  compares  the  differences  in  the 
tangential  stress,  <Jq,  and  a  substantial  reduction  is  noted  in  the 
residual  stress. 

The  usefulness  of  these  results  is  thus  shown.  For  the  quenching 
problem,  the  maximum  quench  time  for  the  desired  metallurgical  phase 
structure  to  be  formed  is  of  interest  because  it  implies  that  the  mate¬ 
rial  will  be  subjected  to  slower  transient  temperatures  and  smaller 
residual  stresses.  Thus,  a  better  understanding  of  the  transient 
temperatures  in  the  quench  tube,  the  resulting  residual  stresses,  and 
the  effect  of  the  quenching  process  is  gained. 
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Figure  1.  Residual  Stresses  in  a  Solid  Cylinder 
Due  to  Material  Transformation  (Trans¬ 
formation  Beginning  on  Outside  Diameter 
and  Progressing  Toward  Center), 
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Figure  3.  Residual  Stresses  in  a  Solid  Cylinder 
Due  to  Quenching. 
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Figure  4.  Residual  Stresses  in  a  Hollow  Cylinder 
Due  to  Quenching. 

Boundary  Conditions: 

hp  h2  =  12.2, - 

h1  =  6.1,  h2  =  12.2,  - 
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